Transfer of optical spectral weight in magnetically ordered superconductors 
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We show that, in antiferromagnetic superconductors, the optical spectral weight transferred to 
low frequencies below the superconducting transition temperature originates from energies that can 
be much larger than twice the superconducting gap A. This contrasts to non-magnetic supercon- 
ductors, where the optical spectrum is suppressed only for frequencies below 2A. In particular, we 
demonstrate that the superfluid condensate of the magnetically ordered superconductor is not only 
due to states of the magnetically reconstructed Fermi surface, but is enhanced by transfer of spectral 
weight from the mid infrared peak generated by the spin density wave gap. We apply our results 
to the iron arsenide superconductors, addressing the decrease of the zero-temperature superfluid 
density in the doping regime where magnetism coexists with unconventional superconductivity. 



I. INTRODUCTION 

Optical measurements reveal crucial spectral and elec- 
tromagnetic properties of a superconductor—^—. For ex- 
ample, the London penetration depth A can be deter- 
mined from the low frequency dependence of the imagi- 
nary part a" (w) of the optical conductivity a (u)): 



a" {lo ^ 0) 
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(1) 



where c is the speed of light. In addition, the super- 
conducting gap A can be obtained from the real part of 
the optical conductivity, a' {uj) , since spectral weight is 
transferred from energies below 2A to the w — > contri- 
bution c^\~^5 (oj) /4. Kramers-Kronig transformation of 
this 5- function term then yields Eq.(IT]). These effects can 
be illustrated if one starts from a Drude conductivity 



ct'H 
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in the normal state, with plasma frequency lo'I = 
^Tre^n/m* and scattering time r. Here, m* is the op- 
tical mass and n is the electron density. In the super- 
conducting state, the entire weight of the Drude con- 
ductivity is transferred to the ^-function if A is larger 
than T^^, leading to the result of the BCS theory A"^ ~ 
ijJi/c^ = ATie^n/ {m*c?^ for clean superconductors'*. In 
the dirty limit. At ^ 1, the transfer of the spectral 
weight below uj — 2/S. can be approximated by 2 A x cr' (0) , 
yielding the well known result for dirty superconductors'^ 
A^2 — ^ I^T(Ji / (? . Consequently, the superfluid conden- 
sate 



dTre-^ 



-A' 



(3) 



is reduced compared to the particle density, Us/n — 
4Ar/7r. The determination of A from the optical spec- 
trum is most efficient for At < 1. Furthermore, the 
investigation of the optical conductivity reveals crucial 



information in unconventional superconductors . In the 
cuprate superconductors, the /-sum rule 



w. 



a' (w) dui 
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was used to analyze whether the anomalous redistribu- 
tion of spectral weight below the superconducting transi- 
tion temperature T^ reveals information about the change 
of the kinetic energy, or more precisely of the optical 
mass m* , upon entering the superconducting state^i^. Fi- 
nally, fine structures in the optical spectrum were used 
to determine the mechanism of superconductivity in the 
cuprates^"— . 

In the recently discovered FeAs superconductor a^^i^'^ , 
the interplay of collective magnetic degrees of freedom 
and superconductivity has attracted great interest, in 
particular given the strong evidence for an electronic 
pairing mechanism with s"' -pairing state^^. In this 
state, the superconducting order parameter has opposite 
signs in different sheets of the Fermi surface separated 
by the magnetic ordering vector Q. In distinction, in the 
conventional s^'^-pairing state, that is expected to origi- 
nate from electron-phonon coupling, the superconducting 
order parameter has the same sign everywhere. The re- 
cent observations of the magnetic resonance modei^, of 
the microscopic coexistence between magnetic and super- 
conducting order^^, and, in particular, of the integer and 
half-integer flux-quantum transitions in a niobium-iron 
pnictide loop^^ give strong evidence for s"' -pairing. 

Important insights about the magnetic, superconduct- 
ing and normal states have also been obtained in mea- 
surements of a {uj]^"—. An analysis based on Eq.([T]) 
and on the Ferrell-Glover-Tinkham (FGT) sum rule^i^ 
(see Eg. p^ below) led to results for the penetration 
depthi2ii2^ that are consistent with the values obtained by 
other techniques (see, for example Ref<^). For instance, 
in Bao.6Ko.4Fc2As2, which has Tc = 37K, the authors of 
RefJ^ found A ~ 2000A at T = lOK. A typical value for 
the largest superconducting gap in the same compound 
was estimated^^ as A ~ lOmeV. Other investigations re- 
solved the individual gaps on the various Fermi surface 
sheets^^"— . 




Figure 1: Transfer of optical spectral weight for a magnetic 
superconductor, considering particle-hole symmetry and a fi- 
nite scattering rate t~^ ~ 0.07Aaf. In the pure antiferromag- 
netic state (red), the spectrum has a peak at Eg = 2Aaf and 
no Drude peak. In the coexistence state (yellow), a fraction 
A / (Aa_f -I- a ) of the total spectral weight of the antiferro- 
niagnetic phase is transferred to a finite Drude peak, while an- 
other fraction A%p/ (Aaf -|- A^) stays at energies above the 
new optical gap Eg = 2\/A^p -I- A^. Thus, a system with- 
out free carriers in the normal state acquires zero frequency 
spectral weight through transfer of weight over energies larger 
than 2 A. 



In the parent compounds BaFe2As2 and SrFe2As2, 
which do not show superconductivity, long range antifer- 
romagnetic order below the Neel temperature T/v leads 
to the formation of a mid- infrared (MIR) peak-*^^, with 
spectral weight being transferred to w ~ 2Aaf due to 
the opening of a spin density wave gap. Here Aaf is the 
single-particle gap for momentum states that are Bragg 
scattered by the magnetic ordering vector (see figure [5]) . 
For BaFe2As2, 2 Aaf ^ lOOOcm-i (i.e. Aaf ^ 62meV), 
yielding Aaf/ (ksT^) — 5, which is not unrealistic for 
itinerant antiferromagnets. For T < Tjv, the low energy 
optical response of the parent compounds is characterized 
by the Drude form, Eq.([2]), however with a significantly 
reduced plasma frequency Wp af — (0.1 — 0.2) w^, where 
ujp is the plasma frequency above T^-^^. Finally, sys- 
tematic investigations of a' {to) in Ba(Fei_xCoa;)2As2 as 
function of temperature and carrier concentration were 
performed in Refsi ^^'^^ . Upon doping, the MIR peak 
gets weaker, consistent with the decrease of the ordered 
magnetic moment M with doping^^. Indeed, assuming 
Aaf oc M and using the results from Ref.^^ for the 
doping dependence of M {T — 0), one finds for doping 
concentrations where T/v — Tc that Aaf — A, strongly 
supporting the view that the same electrons that un- 
dergo Cooper pairing form the ordered moment. This 
is consistent with the recent analysisL^- of the phase dia- 
gram of Ba(Fei_2;Co2;)2As2: neutron scattering exper- 
iments showed that magnetism and superconductivity 
compete strongly, to the extent that the staggered mo- 



ment is suppressed22i2^ below Tc and the magnetic phase 
boundary T/v (x) is bent back towards smaller x-values 
for Tjv < Tc- Theoretical arguments then demonstrate 
that this coexistence is only possible for an s^ -pairing 
state. The optical properties in the regime of simultane- 
ous magnetic and superconducting order are promising 
as they might reveal important information about the 
interplay between the supcrfluid condensate, the normal 
state Drude peak and the MIR peak. 

In this paper we analyze the optical conductivity in the 
magnetically ordered phase as well as in the regime of 
simultaneous magnetic and superconducting order. Us- 
ing parameters suitable to the description of the phase 
diagram-^- of Ba(Fei_a;Coa:)2As2, we numerically obtain 
the optical spectrum below T^. Due to the partial gap- 
ping of the Fermi surface in the itinerant antiferromag- 
netic state, we obtain a Drude peak in addition to a MIR 
peak at a; ~ 2 Aaf, in qualitative agreement with the 
experimental data. Thus, our results give further strong 
evidence for the itinerant character of the magnetically 
ordered state. 

As shown by nuclear magnetic resonance and 
muon spin rotation experiments^^—, antiferromag- 
netism and superconductivity coexist homogeneously in 
Ba(Fei_aCoa;)2As2 for the doping range 0.035 < a: < 
0.059. Recent tunnel diode resonator measurements 
showed that, in this regime, the T = superfluid density 
is reduced when compared to its value in the pure su- 
perconducting state^^. One might expect, at first glance, 
that such a reduction is due to the suppressed plasma 
frequency below T/v, and the inverse squared penetra- 
tion depth A^^ is given by the reduced value w^^p/c^. 
In contrast to this expectation, we find that the super- 
fluid condensate of a superconductor with magnetic long 
range order, while reduced compared to the case with- 
out magnetic order, has X~^ values that are significantly 
larger than w^ af/*^^- ^^ particular, we find a sizable con- 



densate fraction even in the limit where w^ 



af 



0. By 



analytically investigating the simple but relevant limit of 
particle-hole symmetry, we demonstrate that, in the mag- 
netically ordered superconducting state, spectral weight 
with energies ut ~ 2Aaf is transferred from the MIR- 
peak into the singular S (w) term of a' (w), enhancing Hs 
(see figure [1} . This spectral weight is transferred from 
regions of the spectrum that can easily be larger than 
2 A, reflecting the fact that the rigidity of the supercon- 
ducting ground state with respect to transverse current 
fluctuations, while smaller than in the paramagnet, is still 
larger than what the low frequency Drude weight would 
suggest. 

Our results show that the superfluid density in a mag- 
netic superconductor is not only related to the remain- 
ing electronic states of the magnetically reconstructed 
Fermi surface at T = 0, but also to the transfer of spec- 
tral weight around the MIR peak. These conclusions 
are consistent with recent theoretical investigations from 
Vorontsov et al.^^ showing that superconductivity is able 
to coexist with magnetism even when the reconstructed 



Fermi surface at T = is completely gapped. 

The paper is organized as follows: in Section II we re- 
view the basic properties of the optical spectrum of clas- 
sic superconductors. In Section III we present our results 
for the optical conductivity in the magnetically ordered 
phase of the iron arsenides. Section IV is devoted to 
the investigation of the optical spectrum of the coexis- 
tence state and its relationship to the supcrfluid density. 
Section V brings our conclusions and in Appendix A we 
present an explicit calculation of the penetration depth 
using an alternative approach. 



II. OPTICAL CONDUCTIVITY IN 
SUPERCONDUCTORS 

Within the Kubo formalism the longitudinal optical 
conductivity is given by 



cr(cj) 



iO^ 



n(q 
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(5) 



with longitudinal current-current correlation function 

n(q,L.) = -*/ e'"*([j„(q,t),j„(-q,0)]_)di. (6) 
Jo 

where ja is the a-th component of the current operator. 
For the real part of the optical conductivity, we have: 



a' {uj) = DS (w) + a' (uj) , 



(7) 



with the regular contribution a^cg i^)-, that does not con- 
tain a 6 (to) contribution. From Eq.Q, it follows that 
the Drude weight D of the optical conductivity is given 

by 



D ^ -{col + iTTlliq^ 0,u: -^ 0)) 
while the regular contribution is 



H = - 



Imn(q = 0,0;) 



(8) 



(9) 



The current-current correlation function also deter- 
mines the London penetration depth via 



A" 



-2/2 



47rn (q ^ 0, w = 0)) (10) 



i.e. we consider the static current response at small but 
finite momentum, in distinction to the weight D that 
measures the homogeneous (q = 0) response at small u. 
In general, the order of the limits q = 0, uj — )■ ver- 
sus oj = 0, q — >■ 0, matters and yields different results. 
Yet, in the case of a system with gapped excitation spec- 
trum, it was shown in Refj^ that the order in which 
these limits are taken is irrelevant. Thus, in the case of 
a fully gapped superconductor, it follows generally that 
the Drude weight in the superconductor 



D 



c 
4A2 



(11) 



is determined by the penetration depth and thus by the 
superfluid condensate 



(12) 



Aire'' 



Formally, a contribution D ^ in Eq. ([7]) is not a proof 
for superconductivity and may occur in a metallic system 
that is unable to relax its momentum. Then, the metal 
becomes a perfect conductor, where charges are freely 
accelerated by an external electric field. However, in a 
realistic system one always expects scattering events that 
allow for momentum relaxation. In case of a perfect con- 
ductor, such events broaden the singular Drude peak e.g. 



S{oj)^-- ^ 
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with scattering time r. Thus, formally, D = and the 
Drude response becomes part of the regular contribution 
to the conductivity (y'^^„ (w). This is different for a super- 
conductor, where scattering events may cause a reduction 
of the value of D in Eq.Q but do not change the 8 (w)- 
form of the zero frequency contribution. The Meissner 
effect requires that iJi > — 47rn (q — )■ 0, a; = 0). Together 
with the fact that the order of limits does not matter for a 
gapped system^i, follows I? > 0, i.e. the 8 (aj)-form is ro- 
bust. This preservation of the singular lo — Q response is 
a consequence of the unique rigidity of the superconduc- 
tor with respect to transverse current fluctuations. For- 
mally, this rigidity of the superconducting ground state is 
reflected by the smallness of|n(q— >■ 0,w = 0)| compared 

to ^2/ (4^)^ 

A quantitative determination of the penetration depth 
in a superconductor can be performed by analyzing the 
spectral weight transfer in o (uS)^ and is expressed by the 
Ferrell-Glover-Tinkham (FGT) sum rule^i^: 



c-^ ./o+ 



(o-L(^)-'^L(w))dw. 



(14) 



This sum rule relates A to the change of total spec- 
tral weight between the normal state (ns) and the super- 
conducting {sc) state for a; > 0. It follows from the /- 
sum rule, Eq. (J4]), and the emergence of the Db (aj)-term 
only below T^- Thus, upon entering the superconducting 
state, spectral weight is transferred from finite frequen- 
cies to the (5- function at oj = 0. We mention that for 
Eq. (fT4|) to hold, one assumes that the expectation value 



of the optical mass to*, i.e. the value of w^ — Aire n/m* , 
is unaffected by the onset of superconductivity. In the 



.18.22 



FeAs superconductors, this seems to be the cas( 
distinction to the evidence for violation of the FGT sum 
rule in cuprate superconductorsS.. 




Figure 2: Schematic representation of the band structure con- 
sidered here (left panel). The Fermi surface has an elliptical 
electron pocket (blue ellipse) displaced by the magnetic or- 
dering vector Q from the circular hole pocket (red circle). In 
the right panel, we present the Fermi surface at zero doping. 
For simplicity, we displaced the centers of the bands to make 
them coincide. 



III. OPTICAL SPECTRUM IN THE 
ITINERANT ANTIFERROMAGNETIC PHASE 

A. Microscopic model for competing magnetic and 
superconducting order 

To develop a microscopic model of the interplay be- 
tween superconductivity and magnetism, we use a few 
basic ingredients to describe the main features of the iron 
arsenidesi^: the electronic structure is characterized by 
two sets of Fermi surface sheets, a circular hole pocket 
around the center of the Brillouin zone and an elliptical 
electron pocket shifted by the magnetic ordering vector 
Q. The non- interacting part Hq of the Hamiltonian is 
then given by 



■Ho = ^ (6,kcLck<T + ^.k^L'^kff 

kcr 



(15) 



We consider only one hole band located in the center of 
the Brillouin zone with dispersion ^i k, and one electron 
band, shifted by Q from the hole band, with dispersion 
6,k (see figure [5]): 



Ci,k — £1,0 - Ti M 

2m 



C2,k-|-Q 



-£2,0 + TT^ + TT^ - fi, 



2mr 2m. 



(16) 
(17) 



V 



A magnetic interaction / and an interband pairing in- 
teraction V lead to the possibility of antiferromagnetic 
order with antiferromagnetic gap 



^AF = :j XI ^ (cLdk+Q'y 



and superconductivity with coupled gap equations 

^2 = -]^E(^ktC-k;, 



(18) 



(19) 



^1 - -^Y^i^^k+Qt'^-k-Ql 

k 



(20) 



For V < 0, as it would be the case for phonon 
mediated superconductivity, one obtains the s"*""^ state 
(A1A2 > 0), whereas for T^ > it follows the unconven- 
tional sign-changing s"' state (A1A2 < 0). Introducing 

the Nambu operator ^k = ( c^^, cl^. , , dk+Qti '^-k-oj, 
we consider the mean field Hamiltonian 



H = 



I <£k*k 



(21) 



with simultaneous antiferromagnetic and superconduct- 
ing order. Here 



£k 



Ci,k Ai Aaf 

Ai -$i,k Aaf 

Aaf ^2,k+Q A2 

Aaf A2 — C2,k+Q 



(22) 



For details of this model, see Ref a^^i^^ . Now, the q = 
current-current correlation function is given by 

liafs (iUJn) = e^T ^ tr (vkaGw {iVm + iUJn) VwpGk (i'^ni)] 

k,rn 

(23) 
with Gj, {iCiJn) = itj„l — £k and the ve- 
locity matrix in Nambu space Wka = 
c'/9/cadiag(^i,k,6,k,6,k+Q,C2,k+Q)- Here, the in- 
dices a and /3 refer to Cartesian components of vectors, 
ujn — 2mTT is a bosonic Matsubara frequency and 
^'„J = {2m + 1)ttT is a fermionic Matsubara frequency. 



B. Magnetically ordered phase without 
superconductivity 

First, we investigate the optical properties of the pure 
antiferromagnetic state. In this case, one has 2x2 ma- 
trices in Nambu space and the Green's function is given 
by: 



Gsk (i^n) = («W„ - Si_k) (ioJn - E2M, 

i^n — C2,k+Q -sAaf 
-sAaf i^n - Cl,k 



(24) 



where s denotes the spin and i?a,k, the quasiparticle en- 
ergy: 




To evaluate the current-current correlation function, 
Eq. (P5)) . we use the Kramers-Kronig relations: 



GsU {ioJn) 






(26) 



In order to obtain more realistic results, we replace one 
of the iO+ convergence factors above by a finite single- 
particle lifetime ir~^. Then, the real part of the conduc- 
tivity is given only by the regular contribution ^. A 
straightforward calculation leads to: 



c^L i^) = 2e vp 



V^ V^ fa {Ea,ki^) 



a=l k 
{2UJ + EaM — Ea,k) 



{Ea,l 



En 



aM) 



where v^? is the hole-band Fermi velocity and: 



(27) 



fa {Ea,k,^) 



np (EaM) ~ ^F {EaM 
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'-^i.k 



+ C. 



(a) 
2.k 



(28) 



.(a) 



with Cl'^^ = {EaM - Cj,k) {EaM + '^ - C»,k), Fermi func- 
tion np and relative electron-band mass ifia = ma/m. 
Here, we introduced the index a defined as a = 2 for 
a ~ \ and a ~ 1 for a — 2. 

In Ref.^^, we introduced the band structure parame- 
ters that provide a good agreement between the model 
of the previous subsection and the neutron diffraction 
data on Ba(Fei_xCoa;)2As2. Specifically, they are given 
by £1,0 = 0.095 eV, £2,0 = 0.125 eV, m = 1.32mcioctron, 
rux = 2m, and ruy = 0.3to, together with the electronic 
interaction / = 0.95 eV and the assumption that each Co 
atom adds one extra electron. In this subsection, since we 
are only interested in the pure antiferromagnetic phase, 
we set V — Q. In figure |31 we use these parameters to ob- 
tain the real part of the optical conductivity at approxi- 
mately zero temperature and for different Co doping con- 
centrations, with no superconductivity involved. Note 
that the only free parameter here is the single-particle 
lifetime r, which was chosen to be r"^ = 5.5 meV for all 
doping concentrations. This is the same order of mag- 
nitude of the scattering rate associated with the narrow 
Drude peak observed in optical experiments^^. 

Our objective here is not to describe all the details 
of the observed optical spectrum, which would require 
a description of all five Fe-3d orbitals (see, for example, 
Ref."^?.), but rather to understand its main features. From 
figure [3J we see that, in general, the optical spectrum in 
the antiferromagnetic phase has a Drude peak as well as 
a finite-frequency peak. The latter is located very close 
to w ~ 2Aaf, and is associated with the opening of the 
spin density wave gap. As shown in the same figure, Aaf 
only partially gaps the Fermi surface, resulting in a finite 
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Figure 3: Real part of the optical conductivity a' (in units 
of (70 = 2-K^^ pe^v%^/ eo) as function of the frequency to (in 
units of nieV) at T ~ for different Co doping concentra- 
tions: (a) s = 0, (b) I = 0.025, and (c) x = 0.057. We use 
the band structure parameters that consistently describe the 
phase diagram of Ba(Fei_iCoa:)2As2 (see main text for more 
details); only magnetic ordering is considered. The dashed 
gray line refers to the frequency 2Aaf and the inset shows 
the magnetically reconstructed Fermi surface (see figure [2]). 



Drude peak proportional to r. Recall that, for nested 
bands and t — )■ 00, the reconstructed Fermi surface is 
completely gapped and the optical conductivity in the 
magnetically ordered phase vanishes for lo < 2Aaf (see 
Eq.ET] below for A = 0). 

The existence of a peak at 2 Aaf, combined with the 
/-sum rule, Eq. (|3]), implies that the plasma frequency 
associated with the remaining Drude peak in the antifer- 
romagnetic state must be smaller than the plasma fre- 



quency of the Drude peak in the paramagnetic phase, as 
seen experinientallji22, (assuming that the optical mass is 
the same in both situations). Note that the theoretical 
value of Aaf in the undoped sample is Aaf ~ 51 meV 
(Fig. [5^) , which is very close to the value extracted from 
the measured optical spectrum^^. As doping increases 
and the magnitude of the gap decreases, the MIR peak 
gets weaker and moves towards lower frequencies (Fig. 
[5)3), until it is almost completely masked by the Drude 
peak (Fig. [3]:). These results are in general agreement 
with optical measurement a^^'^^i^^ on Ba(Fei_a;COa;)2As2, 
demonstrating the itinerant character of the magnetically 
ordered state in these compounds. 



and we obtain for the optical conductivity 



D = 



4 AIf 



(34) 



For the non-superconducting antiferromagnet (A = 
but Aaf 7^ 0) follows D — 0. This is a consequence 
of perfect nesting that leads to a fully gapped antifer- 
romagnetic state. In the non-magnetic superconductor 
(Aaf = but A ^ 0) the current-current correlation 
function vanishes at T = 0. Then, it follows D = i^p/4: 
and Eg. dTT]) yields the BCS result for the penetration 
depth Ao = ujp/c^- In the general case we obtain for the 
penetration depth 



IV. OPTICAL SPECTRUM IN THE 
MAGNETICALLY ORDERED 
SUPERCONDUCTING PHASE 

So far we have considered only the magnetically or- 
dered state. However, for Ba(Fei_a;Coa;)2As2, super- 
conductivity coexists with antiferromagnetism at very 
low temperatures for 0.035 < x < 0.059. In order 
to achieve a more transparent insight about the opti- 
cal conductivity in the coexistence state, we investigate 
analytically the limit of particle-hole symmetry, where 
Eq = £1,0 = £2,0i i^x = 'Tij, = m and fi = 0. In this case, 
^k = '^i.k — — ^2,k-i-Q, implying that the hole and electron 
Fermi surfaces are identical (perfect nesting) . In the case 
of s^ -pairing we have 



A2 



£k = S,kr3a-3 + AafToCTi + Aricr^ 



(29) 



where A — Ai = — A2 and Tq and ap are the Pauli ma- 
trices that act in Nambu and band space, respectively. In 
case of s++-pairing, we replace tr^ by cfq in the last term. 
Eg. (1291) leads to the single particle Green's function 



Gk (iujn) 



lUJnToao + £k 

"c.2 + A2 + Aip + er 



(30) 



At zero temperature the q = current-current corre- 
lation function is: 



2 J (47r) ^ 



(31) 

where vq = tqGz contains the proper sign of the current 
vertex. Performing the trace over the band and Nambu 
degrees of freedom, we find 



87r2 



.e2+A2-Aip-r!(w + r!) 



nH-^/^^^g - - A^^,^^, 



1 



X- 



^2 -I- Aip + A2 + (w + ny 

In the limit cj ^- follows 



n (w ^ 0) 



47r Aip + A2 



(32) 



(33) 



AIf 



A2 



(35) 



We point out that this result is the same as in the 
case of a charge density wave state coexisting with a con- 
ventional s-wave state^^. As shown in Appendix A, we 
obtain the same result for the penetration depth by ex- 
plicitly analyzing Eq. ljlOp . i.e. by first taking uj ~ and 
then q -^ 0. In this context, the existence of a finite 
n (w = 0, q -^ 0) at T = is related to the fact that one 
of the two coherence factors is not identically zero, in con- 
trast to what happens for non-magnetic superconductors. 
Thus, the rigidity of the non-magnetic BCS ground state 
with respect to transverse current fiuctuations is reduced 
in the magnetically ordered state. 

Note that, formally, the coexistence between supercon- 
ductivity and magnetism is only marginal for particle- 
hole symmetry, as we discussed elsewher&i£. Yet, small 
perturbations in both the chemical potential and the 
ellipticity of the electron band are able to place the 
system in the coexistence regim o^^'^'^'^^i'^^ . Following 
Refsi^i^, we can investigate the effect of these small per- 
turbations on our result for the penetration depth ([35)) 
by considering the perturbed band structure ^2.k+Q = 



-Ci,k — 25^, with 5^ 



82 cos 2if such that 5q 



^^+ — y^- 2m^my ) aiid 62 - — [^ „^„^ j . Here, 
is the angle on the elliptical electron pocket. A straight- 
forward calculation leads to: 



^ 



X-2 



AF 



A2 



4Aip 



3(Aip 



A2)- 



(Si) 



(36) 



with ((52) — (5g+(5|/2. Thus, both perturbations in the 
chemical potential and in the ellipticity lead to a decrease 
in the penetration depth, i.e. to an increase in the value 
of D. Therefore, we can interpret the particle-hole sym- 
metric result (|35|) as an "upper-bound" for A. With this 
in mind, even though the band structure of the pnictides 
is not particle-hole symmetric, it is instructive to substi- 
tute in Eq. (135)) the T = values of A and Aaf obtained 
by numerically solving the gap equations with the band 
structure parameters of Ba(Fei_2.Co2;)2As2 (see Section 
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Figure 4: Ratio of the superfiuid condensate na and 
the electron density n with parameters relevant for 
Ba(Fei-a:Coa;)2As2 in the clean limit, as function of x and 
for T = 0. For x < 0.059 simultaneous antiferromagnetic 
and superconducting order sets in while for x > 0.059 no long 
range magnetic order exists. The inset shows the T = values 
of A (dashed red line) and Aaf (solid red line) as function of 
x, according to the theory of Ref.— . For simplicity we con- 
sider the particle-hole symmetry expression (|35|) and use the 
value of the SC gap referent to the electron band. 



Ill) . The results are displayed in figure U) Remarkably, 
similar values for the relative increase of At=o have been 
recently measured by Gordon et al^ along the coexis- 
tence region of Ba(Fei_a:Coa:)2As2 using the tunnel diode 
resonator technique. 

Next, we analyze the optical conductivity at finite fre- 
quencies. To this end, we evaluate the integrations in 
Eq. (|32p and perform the analytical continuation to the 
real frequency axis, iuj ^ uj + iO"*", obtaining: 



reg 



(c.) 



= < w: 







yf^^^ 



UJ < Eg 

io> E„ 



with the optical gap 



E„ 



2A/Aip + A2 



(37) 



(38) 



In the normal state, A = and the optical conductivity 
is nonzero only for w > 2 Aaf- Entering the supercon- 
ducting state, it follows for A < Aaf that there is no 
spectral weight in the normal state for w < 2A. Thus, 
the finite penetration depth obtained in Eq. (|35|) must be 
due to the transfer of spectral weight that involves en- 
ergies above 2 A. Indeed, analyzing the remaining high 
frequency spectral weight, we find from Eq. p7p that 



'reg 



(cj) duj = 



A2 

^AF 



AIf 



■A2 



(39) 



where the factor of 2 accounts for negative frequencies. 
Thus, the total weight of the non-superconducting an- 
tiferromagnet splits in two parts with ratio (A/ Aaf) • 




4 KJA 



Figure 5: Ratio of the superfiuid density Us and the electron 
density n as function of the antiferromagnetic gap Aaf (in 
units of the superconducting gap A) for s"'"^-pairing (dashed 
line) and s^~-pairing (solid line) at T = 0. 



Below Tc a fraction A^/ (A^p -I- A^) is transferred to 
cj = to yield Meissner screening and a finite penetra- 
tion depth. In addition, the fraction A\p/ (A^p -f A^) 
remains at energies above the optical gap Eg, as illus- 
trated in figure [T] Here we use the fact that Eg of the 
non-superconducting antiferromagnetic state and of the 
state with A and Aaf finite is essentially the same due 
to the reduced ordered moment below Tc, see Ref,^^. 

Analogously, one can rationalize this result using 
Ferrell-Glover-Tinkham (FGT) sum rule, Eq. ^. By 
calculating the difference of spectral weight between the 
non-superconducting state and the superconducting state 
with the same value of Aaf, we obtain exactly the Drude 
weight D: 



/•OO 

2 / [a',,g (w, Aaf, A = 0) - a'^^g [oj, Aaf, A)] dco 

(40) 



4 Aip + A2 



Our analysis for particle-hole symmetry clearly shows 
that, in a magnetic superconductor, the transfer of spec- 
tral weight to the Drude peak below Tc < T/v is not from 
u) < 2 A, but from higher frequencies uj < 2-\/A^p + A^. 
Therefore, analyzing our results for the optical conduc- 
tivity of the pure magnetic phase (figure[3]), as well as the 
experimental optical spectrum, the superfiuid condensate 
is formed not only by the reduced remaining Drude peak 
but also by the significant portion of spectral weight asso- 
ciated with the MIR peak. More importantly, even in the 
absence of a remaining Drude peak in the pure magnetic 
state at T = 0, spectral weight can be transferred to a 
d (a;)-term below Tc, yielding a finite superfiuid density. 

Let us briefiy discuss the situation for S++- 
superconductivity and particle-hole symmetric bands. 



Similar calculations lead to the T = penetration depth 



2AA 



A- 



An 



AF 



(Aip 



A2)ln 



Aaf-A 



='+A 



4AA 



AF 



(41) 



Correspondingly, the high frequency contribution has 
the total weight 1 — (Aq/A) . In figure [S] we plot the 
superfluid density Ug (A, Aaf) as function of Aaf/A for 
S++ and s^ -pairing. We see that for Aaf < A holds 
that Us is larger for s+^-pairing compared to s^ , while 
the opposite is true for Aaf > A. Qualitatively, the two 
behaviors are not very different and would not provide 
a sharp criterion to identify the symmetry of the pairing 
state in the iron arsenides. Yet, previous analysi o^^'^'^i'^^ 
demonstrated that the conventional s^^ state is unable 
to coexist with itinerant magnetism in the pnictides. 



V. CONCLUSIONS 

In conclusion, we analyzed the optical conductivity of 
both an itinerant antiferromagnetic state and a mag- 
netically ordered superconductor. For the pure mag- 
netic phase, using the parameters associated to the phase 
diagram^^'^^ of Ba(Fei_2;Coa;)2As2, we were able to iden- 
tify the main features observed in the experimental opti- 
cal spectrurriF^i^^i^. In particular, for the undoped com- 
pound, we found a reduced Drude peak associated to 
the remaining reconstructed parts of the Fermi surface, 
as well as a mid-infrared (MIR) peak at oj ~ 2Aaf ~ 51 
meV, associated to the gap opened at momentum ko that 
are Bragg scattered by the magnetic ordering vector Q, 
i-e. ^i,ko = C2,ko+Q- Upon doping, the spin density wave 
gap is reduced and, consequently, the MIR peak becomes 
weaker and more masked by the Drude peak. 

The experimentally observed optical conductivity has 
other particular features that are not contemplated by 
our two-band based model, such as high-frequency inter- 
band transitions, a seemingly doping-independent inco- 
herent contribution with a rather long tail, other possi- 
ble low-weight Drude-like peaks at finite frequencies^!^ 
and, of course, the origin of the scattering processes that 
lead to a finite lifetime r. Clearly, a detailed description 
of the optical spectrum has to take into account the ef- 
fects of the other bands that do not participate in the spin 
density wave state"^^, and possibly the role played by dif- 
ferent orbitals that cross the Fermi level. Yet, our simpli- 
fied model that provides a very satisfactory descriptioriF^ 
of the phase diagram of Ba(Fei_2,Coa;)2As2, is able to 
correctly capture not only the main qualitative features 
of the spectrum, but also the order of magnitude of the 
frequency associated to the MIR peak. 

Most interestingly, our results clarify how spectral 
weight is transferred in magnetic superconductors below 



the superconducting transition temperature, even in case 
where the Fermi surface of the ordered magnet is fully 
gapped. In classical superconductors, only the spectral 
weight below the optical gap Eg — 2A is transferred to 
the Drude peak. However, in the case where itinerant 
magnetism is also present and particle-hole symmetry 
holds, the optical gap is given by Eg = 2a/A^p -|- A^, 
involving energies that can potentially be much larger 
than the superconducting gap. Thus, in the regime where 
magnetism and superconductivity coexist in the iron ar- 
senides, the remaining Drude peak of the antiferromag- 
netic phase, whose plasma frequency can be significantly 
smaller than the plasma frequency of the paramagnetic 
state, is not the only origin for the value of the super- 
fluid condensate. Instead, spectral weight associated to 
the higher-frequency MIR peak is transferred to the 5- 
function at w = 0, enhancing the superfluid conden- 
sate. Yet, this superfluid density is always smaller than 
its value in the non-magnetic superconducting phase, in 
agreement with experiments^!. 

This transfer of optical spectral weight is a consequence 
of the unique rigidity of the superconductor with respect 
to transverse current fluctuations. It implies that, even in 
the limit where the pure T = antiferromagnetic phase 
has no Drude peak, it is still possible to obtain a flnite 
superfluid density below Tc- The superfluid density of 
the coexistence state is not only associated to electronic 
states from the remaining Fermi surface, what allows the 
superconducting transition to take place even when a 
Fermi surface would not be present in the magnetically 
ordered state. This non-trivial observation is corrobo- 
rated by recent calculations of Vorontsov et al^, that 
found coexisting itinerant magnetism and s^ supercon- 
ductivity in cases where the reconstructed Fermi surface 
would be completely gapped at T = 0. As they pointed 
out, in these situations both the antiferromagnetic and 
superconducting phases are "effectively attractive" and 
cooperate to form the coexistence state. 
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ful discussions and for sharing their penetration depth 
data prior to publication. This research was supported 
by the Ames Laboratory, operated for the U.S. Depart- 
ment of Energy by Iowa State University under Contract 
No. DE-AC02-07CH11358. 



Appendix A: Calculation of the penetration depth in 
the coexistence region 

Using Kramers-Kronig relations, the current-current 
correlation function (j23l) at flnite momentum and flnite 
frequency can be written as: 



f^ diy r dv' L ImGk+qji^^+^^ ImGi,+q(z. + zO+) 
IIai3 (q, JWn) = e T y / — / tr v^+qa — : — -J—Vkp : 

It is straightforward to evaluate the Matsubara sum. Setting a;„ = yields: 



(Al) 



n„^(q,a. = 0) = e^^ / -/ — 



Up (v) - Up {v') 



tr f Vk+qalmGk+q (i^' + iO+) Wk/ijImGk+q {i' + iO+) 



(A2) 
where np is the Fermi function. The imaginary part of the Green's function can be calculated directly from Eq. pop: 



Gk (ly + iO^ 



VTQCro + £k 

2^k 



1 



1 



ly- Ek + iO+ iy + Ek + iO+ 



(A3) 



where we defined the positive excitation energy Ey^ — \/^'^ + A^p + A^. It follows that: 



ImGk (^ + *0+) = _!L(^^I^±M [S {oj -E^)-S{u + E^) 



(A4) 



Substituting Eq. (jA4l) in the expression (|A2p . we can evaluate the frequency integrals as well as the trace in Nambu 
space. In the limit of small momentum, we obtain: 



Ha/j (q ^ 0, w = 0) = Vmv^2e^v^vl 



q->0' 



np (-Bk) - np (£^k+q) 



-Bk — i?k+q 



2 a2 

AF 



^ ^ aa+q + A"-A 

-Bk-Ek+q 



np (£'k) - np (-^k+q) 



yielding, for a two-dimensional isotropic superconductor: 



-, _ CkCk+q + A - A^p 
-Bk£^k+q 



(A5) 



na/3 (q ^ 0, W = 0) = -Vpe^Safj ^ 



2 - 



dnp\ 



8E-JV 



^-^,Aip^^^^^r/3i.k 



i?, 



^^ 



V 2 



(A6) 



We introduce the density of states p and take T — 0, obtaining: 



Ha^ (q ^- 0, w = 0) = -pvj.e'^SaiB 



d^ 



25 



^2 + A2 + A2 



AF 



A2 
^AF 



^AF 



A2 
^AF 



ie + A2 



A2 

^AF 



)3/2 

_iA7l 



For A 7^ or Aaf t^ the first term vanishes, whereas 
the second one gives: 



Tlafj (q -^ 0,w = 0) = -2pvle'^Sa,f, 



A2 
^AF 



A2 

^AF 



(A8) 

which leads to the same result as Eq. ((35)) from the main 
text. Notice, from Eq. (IA5I) . that the non-zero value as- 
sumed by the current-current correlation function at T = 

1 — ^^^'^+^ AE . 

^k'^k + q J 

For a non-magnetic superconductor, this term goes to 



zero as q — > for any temperature; then, the only contri- 
bution to the current-current correlation function comes 



Ckik 



hA^-A^, 



from the usual coherence factor ( 1 -, gr^;^ 

\ ^ki!'k+q 

whose prefactor vanishes at T = due to the existence 
of a gap in the quasiparticle energy spectrum. In both co- 
herence factors, the relative minus sign between A2 and 
A^p is a result of the fact that while ^ and A change from 
one Fermi surface sheet to the other, Aaf stays the same. 
Therefore, the change in the penetration depth cannot be 
attributed to a change only in the density of states, in 
accordance to our analysis of the finite frequency optical 



spectrum. 
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